Ekner et al. Parasites & Vectors 201 1, 4:182 
http://www.parasitesandvectors.eom/content/4/1/182 



e Parasites 
&Vectors 



RESEARCH Open Access 



Anaplasmataceae and Borrelia burgdorferi sensu 
lato in the sand lizard Lacerta agilis and 
co-infection of these bacteria in hosted Ixodes 
ricinus ticks 

Anna Ekner 1 , Krzysztof Dudek 1 , Zofia Sajkowska 1 , Viktoria Majlathova 2 , Igor Majlath 2 and Piotr Tryjanowski 3 * 



Abstract 

Background: Anaplasmataceae and Borrelia burgdorferi s.l. are important tick-borne bacteria maintained in nature 
by transmission between ticks and vertebrate hosts. However, the potential role of lizards as hosts has not been 
sufficiently studied. 

Results: The current study showed that 23 of 171 examined sand lizards Lacerta agilis were PCR positive for 
Anaplasmataceae. The nucleotide sequences of the several selected PCR products showed 100% homology with 
Anaplasma spp. found in Ixodes ricinus collected in Tunisia and Morocco (AY672415 - AY672420). 1.2% of lizard 
collar scale samples were PCR positive for B. lusitaniae. In addition, 12 of 290 examined /. ricinus were PCR positive 
for B. burgdorferi s.l. and 82 were PCR positive for Anaplasmatacea. The number of ticks per lizard and the number 
of ticks PCR positive for both microorganisms per lizard were strongly correlated. Moreover, we found a significant 
correlation between numbers of ticks infected with Anaplasmataceae and with B. burgdorferi s.l. living on the same 
lizard. However, there was no significant correlation between detection of both bacteria in the same tick. 

Conclusions: To the best of our knowledge, this is the first report of Anaplasmataceae DNA and additionally the 
second report of B. burgdorferi s.l DNA detection in the sand lizard. 
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Background 

One of the most widespread bacterium transmitted by 
ticks is Borrelia burgdorferi s. 1., an agent of Lyme bor- 
reliosis [1,2]. Reservoirs of B. burgdorferi s.l. are verte- 
brates and special associations between Borrelia strains 
and particular groups of vertebrate hosts have been 
reported [3]. B. lusitaniae was the most common strain 
detected in lizard species and in ticks feeding on them 
[4-6]. 

Ticks are the main vector of other microorganisms, 
such as intracellular bacteria from the family Anaplas- 
mataceae [7,8], which attract the attention of public 
health professionals worldwide. One of the most 
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important species of this family is Anaplasma phagocy- 
tophilum which causes human anaplasmosis (HA), for- 
merly known as a human granulocytic ehrlichiosis 
(HGE) [9,10]. A. phagocytophilum is an obligate intracel- 
lular bacterium infecting the neutrophils of various 
mammalian species [9]. Anaplasmataceae are main- 
tained in nature by transmission between vectors and 
reservoirs. Reservoirs of the bacteria are vertebrates, 
mainly rodents and ruminants [11,12]. The potential 
role of reptiles as hosts or reservoirs is not known and 
has not been sufficiently evaluated. To date, despite 
being found in ticks feeding on reptiles [13-16], Ana- 
plasma spp. has been detected only in three lizard spe- 
cies, Sceloporus ocidentalis, S. graciosus and Elgaria 
coeruleus, living in North America [13]. 

The most common species of ticks in Europe, Ixodes 
ricinus [17] feeds on a wide variety of vertebrate hosts, 
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such as mammals and birds [8,18-22]. It also parasitises 
reptile species, and larvae and nymphs often feed on 
lizards [5,23-26]. /. ricinus may be infected simulta- 
neously with different combinations of bacteria 
[10,27-29]. Observed microorganism co-infection rates 
suggest that the risk of infection with one tick-borne 
bacterium is not independent of other bacteria [30]. 
There are only a few studies on co-infection of bacteria 
in wild vector or host populations [31,32]. In addition, 
as far as we are aware, there has been only one study on 
co-infection of Anaplasmataceae and B. burgdorferi s.l. 
in ticks feeding on lizard species, moreover, on Lacerta 
viridis [16]. 

In this paper we show the role of sand lizard (Lacerta 
agilis) in the transmission cycle of important tick-borne 
pathogens, Anaplasmataceae and Borrelia burgdorferi s. 
1. Moreover, we found a significant correlation between 
numbers of ticks infected with Anaplasmataceae and 
with B. burgdorferi s.l. living on the same lizard. 

Results 

290 ticks (176 nymphs and 114 larvae) were found in 69 
of 171 (40.4%, 95% CL: 32.9 - 48.1) examined lizards. 
All of them were identified as /. ricinus. B. burgdorferi s. 

I. DNA was detected in 12 of 290 ticks (4.1%, 95% CL: 
2.2 - 7.1; Table 1) taken from 10 lizards (5.9%, 95% CL: 
2.8 - 10.5; Table 2). B. lusitaniae DNA was detected in 
8 ticks (66.7%, 95% CL: 34.9 - 90.1), B. burgdorferi s.s. 
DNA in 2 ticks (16.7%, 95% CL: 02.1 - 48.4), the 
remaining two strains were not determined in RFLP 
analyses (Table 1). Among 171 examined lizards (43 
adult females, 59 adult males, 26 sub-adults, 43 juve- 
niles), 2 individuals (1 female, 1 sub-adult; 1.2%, 95% 
CL: 0.1 - 4.1) were PCR positive for B. lusitaniae. Both 
of them had ticks at the time of collection (Table 3). No 
other Borrelia strains DNA were detected in lizard col- 
lar scales. The number of ticks per lizard was strongly 
positively correlated with the number of ticks PCR posi- 
tive for B. burgdorferi s.l. (r = 0.57, n = 171, p < 0.0001; 
Figure 1). 

Anaplasmataceae DNA was detected in 82 (52 
nymphs, 30 larvae) of 290 ticks (28.3%, 95% CL: 23.2 - 
33.8; Table 1) taken from 29 lizards (17.0%, 95% CL: 

II. 7 - 23.4; Table 2). In the body scales of 171 examined 
lizards Anaplasmataceae DNA was detected in 23 



individuals (13.4%, 95% CL = 8.7 - 19.5). Obtained 
sequences of several PCR products showed 100% 
homology to each other. The sequence was compared 
with GenBank entries by Blast N2.2.13 and revealed 
100% homology with Anaplasma spp. strains found in /. 
ricinus ticks collected on vegetation in Tunisia and 
Morocco (AY672415 - AY672420). Among the 23 
lizards (14 males, 8 females, 1 sub-adult) PCR positive 
for Anaplasmataceae in the skin, 15 had no ticks 
(65.2%, 95% CL: 42.7 - 83.6) at the time of collection 
(Table 3). However, lizards that were PCR positive for 
Anaplasmataceae had a higher average number of ticks 
(mean ± SE; 5.0 ± 2.2) compared to lizards PCR nega- 
tive for Anaplasmataceae (1.2 ± 0.2; U-test, Z = -2.51, p 
= 0.012), as well as a higher number of ticks PCR posi- 
tive for Anaplasmataceae (2.2 ± 1.0 vs. 0.2 ± 0.1; U-test, 
Z = -2.68, p = 0.007). The number of ticks per lizard 
and the number of ticks PCR positive for Anaplasmata- 
ceae per lizard were strongly correlated (r = 0.73, n = 
171, p < 0.001; Figure 2). 

Kendall's tau coefficient shows that the number of 
ticks PCR positive for Anaplasmataceae were strongly 
correlated with the number of ticks PCR positive for B. 
burgdorferi s.l. (r = 0.42, n = 171, p < 0.0001) feeding on 
the same lizard. 6 (2.1%, 95% CL = 0.8 - 04.5) ticks con- 
tained both B. burgdorferi s.l. and Anaplasmataceae 
DNA. However, Kendall's tau coefficient showed no sig- 
nificant correlation between detection of B. burgdorferi 
s.l. and Anaplasmataceae DNA (r = 0.10, n = 290, p = 
0.088) in the same tick. 

The co-infection index (Ic) for interactions between B. 
burgdorferi s.l. and Anaplasmataceae was 0.75. However, 
the difference between the number of obtained and 
expected co-infections was not significant {j 2 = 0.02, p 
= 0.89). This suggested that there was no significant 
association between the bacteria. 

We could not analyse a co-infection in lizards, because 
of too the small sample of lizards PCR positive for Bor- 
relia burgdorferi s.l. 

Discussion 

Lizards as hosts of ticks are exposed to various tick- 
borne pathogens. Previous studies have showed that rep- 
tiles are included in transmission cycles of B. burgdorferi 
s. 1. [4,5,33]. In our study, B. lusitaniae DNA was 



Table 1 Number of ticks PCR positive for B. burgdorferi s.l. and Anaplasmataceae. 





ticks collected from the lizards 




290 


ticks infected with Borrelia 


ticks infected with Anaplasmataceae 


12 


82 


ticks infected with 8. lusitanieae ticks infected with 8. burgdorferi s.s. 


ticks infected with 8. burgdorferi s.l. 


8 2 


4 



Ekner et al. Parasites & Vectors 201 1, 4:182 
http://www.parasitesandvectors.eom/content/4/1/182 



Page 3 of 7 



Table 2 Number of lizards with and without ticks PCR positive for Borrelia and Anaplasmataceae. 



lizards without ticks 




lizards with ticks 




102 




69 






lizards infected with Anaplasmataceae 


lizards infected with Borrelia 


lizards non-infected 




29 


10 


30 



detected in 1.2% of lizards. Compared with other papers, 
this is a low infestation [4,5]. However, this is only the 
second detection of the bacterium in sand lizards [4] 
and the third detection of the B. lusitaniae strain in 
Poland [6,34]. In all lizards, as well in most of the ticks, 
B. lusitaniae was detected, which confirms the connec- 
tion of this strain with reptile species [3,35,36]. 

During the study 28.3% of ticks feeding on lizards 
were infected with Anaplasmataceae. This is high com- 
pared to other studies on ticks from lizard species 
[13,14,16]. However, despite Anaplasma spp. DNA 
being previously detected in ticks collected from reptiles 
[14], it had only been detected in three lizard species, 5. 
occidentalis, S. graciosus and Elgaria coeruleus [13]. In 
our study, 13.4% of lizards were PCR positive for Ana- 
plasmataceae. To the best of our knowledge this is the 
first detection of the bacterium in sand lizard, moreover, 
the first among lizards living in Europe, and only the 
fourth among lizards worldwide. In addition, the pre- 
vious report concerned 10.2% of lizards living in Califor- 
nia which were infected with A. phagocytophilum [13]. 
In presented study, some lizards were PCR positive for 
different species of Anaplasmataceae than in the pre- 
vious study on reptiles [13], namely Anaplasma spp. 
100% homology with strains found in /. ricinus collected 
on vegetation in Tunisia and Morocco (AY672415 - 
AY672420) [37]. The pathogen could be transferred to 
Poland together with exotic reptiles or other animals 
[14]. On the other hand, lizards could not be analysed 
for presence of that species before. 

Lizards PCR positive for Anaplasmataceae had more 
ticks than non infected individuals. Moreover, the num- 
ber of ticks feeding on a lizard was strongly correlated 
with both the number of ticks PCR positive for Anaplas- 
mataceae and PCR positive for B. burgdorferi s.l. It may 
result from that the more the lizard has ticks, the more 



Table 3 Number of lizards PCR positive for B. burgdorferi 
s.l. and Anaplasmataceae and presence of ticks in lizards 
with the bacteria during the time of collection. 



lizards cached during the study 
171 



lizards infected with Borrelia 
2 



lizards infected with 
Anaplasmataceae 
23 



chance it has ticks with bacteria. Such correlations sug- 
gest that the more ticks feed on a lizard, the greater 
chance of contact with infected ticks, and hence more 
chance of acquiring tick-borne infection. The results 
confirmed previous studies [38,39] which show that the 
probability of host exposure to a tick-borne pathogen is 
correlated with tick abundance. Therefore, the best way 
to avoid infection with tick-borne disease is to avoid 
areas with a high density of its vectors [18]. 

Ticks can be infected with two or more microorgan- 
isms simultaneously [16,38,40,41], but relationships 
between them in ticks can be varied. Some of them dis- 
play antagonistic interactions, others positive, and many 
evidently do not interact [38]. Interpretation of the 
results can be difficult, because many factors, other than 
simply interaction between microorganisms, may influ- 
ence the number of co-infections. In addition, the 
occurrence of tick-borne pathogens in nature may be 
influenced by a number of factors, such as microclimate 
conditions, vegetation, and tick density [42]. In the pre- 
sent study we did not show any correlation between 
detection of B. burgdorferi s.l. and Anaplasmataceae, 
occurring in the same ticks. The results may suggest a 
lack of interaction between the bacteria. These results 
are in contrast to those obtained by Vaclav [16], where 
Anaplasmataceae had a negative influence on B. 




lizards with ticks 
2 



lizards without 
ticks 
0 



lizards with ticks lizards without ticks 
8 15 



Figure 1 Correlation between total number of ticks and the 
ticks PCR positive for B. burgdorferi s.l.. Positive correlation 
between the number of the ticks feeding on a lizard and the 
number of the ticks PCR positive for B. burgdorferi s.l. (r = 0.57, n = 
171, p < 0.0001). 
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Figure 2 Correlation between total number of ticks and the 
ticks PCR positive for Anaplasmataceae Positive correlation 
between the number of the ticks feeding on a lizard and the 
number of the ticks PCR positive for Anaplasmataceae (r = 0.73, n = 
171, p < 0.001). 



lusitaniae. However, in the same study, B. lusitaniae 
positively influenced Anaplasmataceae prevalence, and 
co-infection of both bacteria in ticks was higher than 
expected [16]. To the best of our knowledge, the present 
work is the first study of mixed infection of B. burgdor- 
feri s.l. and Anaplasmataceae DNA in ticks feeding on 
L. agilis [16]. Knowledge of the multiple infection is 
very important for public health, especially for a correct 
diagnosis and prophylaxis of tick-borne diseases, as well 
as prognosis of mixed infection in humans. Moreover, 
hosts infected by several different pathogens can have 
different symptoms of a disease [43]. Knowledge about 
the local occurrence of pathogens may be useful when 
disease symptoms of patients bitten in that locality are 
unclear [10]. It is important to know how the bacteria 
can coexist in individual ticks as a prerequisite for the 
occurrence of co-transmission from tick to the host 
[30]. 

In our study, 65.2% of 23 lizards PCR positive for 
Anaplasmataceae did not have any ticks at the time of 
collection, which may suggest that this bacterium is 
maintained in a lizard body longer than the source of 
the infection. However, we still do not know if lizards 
are reservoirs of Anaplasmataceae or just have organ- 
isms deposited in them by infected ticks. A previous 
study, where the author experimentally infected lizards 
with the bacterium, concluded that lizards were not 
reservoir hosts for Anaplasma phagocytophilum [13]. 
However, it does not mean that lizards cannot be a 
reservoir of the other species of Anaplasmataceae. 
Moreover, lizards may influence the transmission cycle 



of bacteria in areas where there are significant hosts for 
ticks [2]. 

Conclusions 

To the best of our knowledge, the current study is the 
first report of Anaplasmataceae DNA, and the second 
report of B. burgdorferi s.l DNA detection in a European 
lizard species, namely the sand lizard Lacerta agilis. 
Obtained results suggest that lizards may be a reservoirs 
of this pathogen and can influence the transmission 
cycle of the bacteria in some areas. Moreover, we found 
a significant correlation between numbers of ticks 
infected with Anaplasmataceae and with B. burgdorferi 
s.l. living on the same lizard. This knowledge may be 
important in the estimation of the dispersion of the 
tick-borne pathogen and/or sources of potential human 
infection. 

Materials and methods 

Study area and the study species 

The study was carried out in March - September in 
2008 and 2009 in an extensive farmland area in the Bar- 
ycz valley, in Poland (51°34'N, 17°40'E, elevation 110- 
170 m). This study area is characterised by intensively 
farmed land with a varied mosaic of arable fields, mea- 
dows, small woodlots and scattered trees and shrubs of 
different ages, dominated by white willow Salix fragilis, 
silver birch Betula pendula, black poplar Populus nigra 
and pine Pinus silvestris. It contains both dry sandy 
areas and moist areas (for details see reference [44]). 

The sand lizard is a short-legged, rather robust, small 
to medium sized lizard (up to 110 mm snout to vent 
length (SVL)) from the family Lacertidae. It is a ground- 
dwelling and strongly diurnal species with one of the 
widest distribution ranges of all reptiles [45]. In the 
study area the sand lizard is a common species, and an 
average of 0.37 individuals were noted on 200 m trans- 
ect route [44]. 

Lizard and tick sampling 

Lizards were captured using landing fishnets or by hand, 
then aged (adult, sub-adult and juvenile) and sexed. Ani- 
mals were examined for the presence of ticks, which 
were removed with forceps and stored in 70% ethanol. 
Ticks were identified to species and aged using a bino- 
cular microscope, according to Siuda [17]. 

From each individual lizard a skin biopsy was taken 
from collar scales (3-4 mm in length) with sterile scis- 
sors and put in separate vials with 70% ethanol. This 
method had been previously successfully used to detect 
tick-borne pathogens in reptiles [4,5]. The collar is an 
extension of the skin, hence this method is only mini- 
mally invasive to the lizard. The sample is also easy to 
obtain. To avoid resampling the same individual, lizards 
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were permanently marking using Medical Cautery Units 
(unpublished observations). 

Lizard capture was carried out according to Polish law 
and the ethical commission for the study on animals 
(LKE 12/2007). 

DNA isolation 

Immediately prior to extraction, ticks and tissues were 
dried for 30 min to evaporate the ethanol. Each sample 
was cut with a disposable sterile scalpel. Genomic DNA 
from lizard scales and from ticks was isolated by alkaline 
hydrolysis, according to previous reference [46], with a 
30 min. incubation time. Cut samples were incubated in 
the presence of 100 ul ammonium hydroxide (0,5 mol/1) 
at 100°C for 30 minutes in 1.5 ml tube, followed by 10 
minutes at 100°C with the tube open. Isolated DNA was 
stored at -20°C. 

Polymerase Chain Reaction (PCR) 

PCR amplification was performed in a total of 25 ul 
reaction mixture of a MasterTaq DNA polymerase kit 
(Eppendorf AG, Hamburg, Germany) containing 10.4 ul 
of deionized water, 5.0 ul of 5 x TaqMaster PCR Enhan- 
cer, 2.5 ul of 10 x Taq buffer (with 15 mM Mg 2+ ), 1.5 ul 
of a 25 mM solution of Mg(OAc)2, 0.1 ul of Taq DNA 
polymerase (5 U/ml), 0.5 ul of dNTP-mix (10 mM) (Fer- 
mentas, Vilnius, Lithuania), 1.25 ul of each primer (10 
pmole/ul) (Invitrogen, Paisley, Scotland), and 2.5 ul of 
DNA template. 

In order to verify that DNA had been successfully iso- 
lated from each tick, primers for the fragment of the 
tick's mitochondrial cytochrome b gene (620 bp) were 
used [47]. Verification of successfully isolated DNA 
from lizard scales was carried out using primers for the 
fragment of the vertebrate's 12S rDNA [48]. Seven nega- 
tive samples of ticks and two negative samples of lizards 
were excluded from further analysis. 

Samples with successfully isolated DNA of ticks and 
lizard were examined for the presence of Anaplasmata- 
ceae DNA by amplifying a portion of the region of the 
16S (rrs) rRNA gene of the family Anaplasmataceae 
[37,49]. Samples were also examined for the presence of 
B. burgdorferi s.l. by amplifying a portion of the 5S 
(rrfA)-23S (rrlB) rDNA intergenic spacer [50]. The PCR 
products were electrophoresed on a 1% agarose gel, 
stained with Gold View Nucleic Acid Stain, and visua- 
lized with a UV transilluminator. 

RFLP analysis 

The positive PCR products of the 5S-23S rDNA inter- 
genic spacer regions were further analyzed by RFLP. 
Previously extracted DNA of B. afzelii, B. garinii, B. 
valaisiana, B. burgdorferi s.s. and B. lusitaniae were 
used as positive controls. For each positive sample 13 ul 



of amplified DNA were digested at 65°C overnight in a 
solution containing 5 U of Trul I (300 u/ml) and 1 x 
Buffer R (Fermentas). Electrophoresis was carried out in 
16% polyacrylamide gel at 150 V for 3 h. The gels were 
stained with SYBR Gold nucleic acid gel stain (Molecu- 
lar Probes, Leiden, The Netherlands) for 20 min, and 
bands were visualized with a UV transilluminator. 

All procedures, DNA isolation, PCR, and electrophor- 
esis were performed in separate rooms using different 
pipettes and racks, with separate lab coats and disposa- 
ble gloves worn in each laboratory to prevent carry-over 
contamination and to avoid false-positive results. PCR 
mixture was prepared in a sterile PCR box. All liquid 
handling procedures were performed using disposable 
sterile filter tips. In each DNA isolation and PCR reac- 
tion, a negative control (water) was included. 

DNA sequencing of PCR products 

Three randomly selected PCR products of 16S rDNA of 
the Anaplasmataceae family were sequenced in the 
Laboratory of Biomedical Microbiology and Immunol- 
ogy at the University of Veterinary Medicine and Phar- 
macy in Kosice. Sequencing only a few selected samples 
to exactly verify a bacterium species is an acceptable 
procedure [14]. Prior to the sequencing, PCR products 
were purified using a QIAquick PCR purification kit 
(Qiagen). The complementary strands of each sequenced 
product were manually assembled. 

Index of co-infection (l c ) 

Ginsberg [38] developed an index of co-infection (I c ), 
which quantify the degree of departure of the number 
of mixed infections from independence. This is defined 
as the difference of the number of co-infections from 
the number expected due to chance alone, as a percen- 
tage of the total number of infected ticks in the sample. 

I c = [(O - E)/N] x 100, 

where: O = number of observed coinfections, E = 
expected number of co-infected ticks due to chance 
alone, N = total number of ticks infected by either or 
both microorganisms. 

E = [(a + b)(a + c)]/(a + b + c + d), 
N = a + b + c, 

where: a = number of ticks infected with both bacteria 
(equals O), b = number of ticks infected only with 
microorganism 1, c = number of ticks infected only with 
microorganism 2, and d = number of ticks not infected 
with either microorganism. Ic is positive when the num- 
ber of co-infections is greater than expected, and nega- 
tive when there are fewer co-infections than would be 
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expected due to chance alone. Significance of the index 
was calculated by a chi-square test. 

Statistical analysis 

To improve sample size and show more general pat- 
terns, data from the two breeding seasons (2008 and 
2009) were pooled. Statistics were performed using 
SPSS for Windows, and all tests are two-tailed. Confi- 
dence limits (95% CL) for binary, presence-absence, data 
were calculated in an Excel macro. 
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